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Abstract 

A method is presented for modelling the optical properties of a pho- 
tonic crystal structure mounted on a substrate which is thick enough that 
the light reflected from the back is incoherent with reflections from the 
front. Transmission and reflection spectra are presented for the cases 
where the structure is a multilayer planar microcavity and an opal pho- 
tonic crystal. The contributions from the incoherent reflections can be 
very large, and in the photonic crystal case, change the positions of the 
spectral peaks. 

There are many situations where a photonic crystal structure which is stud- 
ied experimentally is mounted on a thick substrate. In modelling the optical 
properties of the actual structure, wc need to consider the effects of interference 
between light reflected off different interfaces. However, if we also allow inter- 
ference between waves reflected from the structure and off back of the substrate, 
the spectra which arc obtained are very unrealistic, being dominated by closely 
spaced Fabry-Perot fringes. This problem is typically avoided in calculations 
by treating the substrate as semi-infinite, eliminating any reflections. In reality, 
the interference effects arc not seen because inhomogencities in the substrate, 
such as thickness variations, mean that any waves which has passed through the 
substrate is effectively incoherent with the front reflections. In this letter, we 
present a simple but much more realistic theoretical treatment which incorpo- 
rates multiple reflections in the substrate as incoherent waves. Our results show 
that these reflections can make very significant contributions to the spectra, so 
the semi-infinite substrate approach is not appropriate for accurate comparisons 
with experimental data. 

The method we describe is very adaptable, and can be applied to any struc- 
ture for which the intensity transmittance and reflectance on a semi-infinite 
substrate can be calculated, either analytically or numerically. As indicated in 
Fig.(l), we call these quantities Tf and i?/ for light passing from air into the 



Figure 1: Definitions of tlie front and back transmittance and reflectance coef- 
ficients. 



substrate, T/ and Rf for the reverse direction. We also require the the trans- 
mittance and reflectance of the back of the substrate, T;, and Rb, which come 
from the standard Fresnel expressions. 

The total transmittance, T, is obtained by summing over an infinite number 
of incoherent reflections from the surfaces of the substrate. If the absorption 
coefficient of the substrate, a, is non-zero there is an attenuation exp {—al) on 
every pass. The result is 



T = T, 



l + RbRfC 



-2al 



(RbRf e 



-2al\2 



The 



= Tf{l-e-^"'RbRf> ±b 
Similarly the total reflectance, TZ, is 

n = Rf + Tf\l + RbRf 

+iRbRf e-^"^f + ■ ■ ■] Rbffe 



(1) 



-2al 



R 



f 



Tj (1 - e-'^'^^RbRf)-^Rbff e 



-2al 



(2) 



Since there is no interference, these results are independent of the thickness of 
the substrate, I, except for the attenuation due to absorption. 

For a transparent substrate, with no structure on the surface, the reflectances 
are all the same, Rf = Rf = Rb = R, and also Tf = Tf = Tb = I — R. Eqs.(l) 
and (2) then reproduce the well known results for incoherent reflections. 
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A further check is the case of a thin fllm on the surface of the substrate. The 
transmittance of such a fllm with incoherent substrate reflections is given by 
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Figure 2: Solid lines: normal incidence transmission (a) and reflection (b) 
spectra for the microcavity structure on a substrate. Dashed lines: results 
obtained when the incoherent multiple reflections in the substrate are neglected. 

SwanepoelfTI. It is straight forward to obtain expressions for Tf and Rf for a 
simple film, and Eq.(l) then correctly reproduces this result 

Turning to the case of a one-dimensional photonic crystal, we next calculate 
the properties of a planar microcavity on a thick substrate. The parameters 
used are typical for a GaAs/AlGaAs structure, with a design wavelength of 
Ao = 900nm. The cavity is a A/2 thickness [3 layer of refractive index 3.5. 
with the mirrors on either side comprising 10 bilayer Bragg stacks with layer 
thicknesses A/4 and high and low indicies of 3.5 and 3. The substrate index is 
3.5 and its absorption coefficient a = 0. 

In this case, it is easiest to obtain the front coefficients, Rf, Rj, Tf and Tf 
from a transfer matrix calculation!?]. We have done this numerically, and the 
transmission and reflection spectra calculated using Eqs.(l) and (2) are shown in 
Fig. (2). The dashed lines on the figure show the results obtained when multiple 
reflections in the substrate are ignored, by putting i?f, — 0; for the reflectance, 
this is simply Rf, the structure reflectance, while for the transmittance it is 
TfTf,. The incoherent reflections always add to T and TZ, and their contribu- 
tion is significant, except where the reflectance of the structure is high, and no 
light passes through the microcavity into the substrate. Where the structure 
reflectance is zero or very small, all the light passes through in both directions, 
and we just see the reflectance and transmittance of the back of the substrate, 
which are ~ 0.31 and 0.69 respectively. 

The situation becomes more complicated when we consider a two or three 
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Figure 3: Normal incidence, unpolarised transmission (a) and reflection (b) 
spectra for the 3-layer opal structure on a glass substrate. The solid and dashed 
lines are results with, and without, incoherent reflections, as in Fig. (2). 

dimensional photonic crystal fabricated on a substrate, so the structure has an 
in-plane periodicity. The effect of this periodicity is to give rise to diffraction, 
so for some wavelengths a single incident wave will produce multiple diffracted 
waves in the substrate, or in both substrate and air. Despite this, Eqs.(l) and 
(2) still apply, though the symbols needs to be reinterpreted as matrices. 

Consider a wavelength at which there are Na waves in air, and Ng waves 
in the substrate, including both straight through and diffracted directions. In 
general, photonic crystals cause mixing of polarisations, so both TE and TM 
waves need to be included and the N values will be twice the number of direc- 
tions. We can then describe the front reflection and transmission processes by 
matrices; for example, Tf becomes a x Na matrix, whose element gives 
the intensity generated in substrate wave i for illumination in exterior wave j. 
We assume again that there is some method available to obtain these matrices; 
in the example presented in Fig. (3), they are calculated using the scattering 
matrix method of Ref. 5 . The back coefficients Tf, and Rt become diagonal ma- 
trices, with elements given by the standard Fresnel expressions for finite angle 
transmittance and reflectance. Knowing these matrix elements, we can add up 
intensities, as in the derivation of Eqs.(l) and (2), with the complication that 
we now need to sum over every possible path, comprising the combinations of 
the different directions and polarisations that can occur on each pass through 
the substrate. This summation is obtained by replacing expressions such as 
RbRf by the equivalent matrix products, since the sum over the internal index 
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in the product replicates the sum over paths. Eqs.(l) and (2) then give the total 
transmittance and reflectance T and TZ, now as Na x A^^ matrices. 

In Fig. (3), we show the results of this procedure for an opal structure on a 
glass substrate. The opal consists of three licp-ordercd layers of 500nm spheres, 
with refractive index 1.5; the substrate also has index 1.5 and absorption a = 0. 
For wavelengths longer than 649. 5nm, there is no diffraction, and the physics is 
similar to the planar microcavity, though we have 2x2 matrices because of the 
polarisation mixing induced by the photonic crystal. The structure reflectance 
is very low, so the substrate contributes just one back reflection, which adds a 
constant 0.04 to the reflectance but does not affect the transmittance. 

At shorter wavelengths, there are six diffracted directions in the substrate, 
but, in the range of the figure, none in air, so A^a = 2 and Ns = 14. Again, the 
structure reflectance is quite low, but the diffracted waves, unable to propagate 
in air, experience total internal reflection at the back of the substrate. They 
then interact with the opal again and are partially re-diffracted back into the 
normal direction, to provide strong contributions to both the reflectance and the 
transmittance. It is notable in the reflectance that these contributions not only 
are much larger than those from the front reflection, but they are also different 
spectrally, with peaks appearing at different wavelengths. It is clear from this 
that it may be very misleading to neglect the substrate reflections, though we 
believe that this has been done in all previous published calculations. 

In conclusion we have shown that, if transmission and reflection spectra 
can be calculated for an optical structure, there is a very simple formulism by 
which to add the contributions due to multiple incoherent reflections occuring 
when the structure is mounted on a finite substrate. We have used this method 
to calculate the optical properties of a planar microcavity structure and an 
opal photonic crystal, and shown that the substrate reflections can modify very 
significantly the spectra which are obtained. 
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Supplementary Information: Derivation of trans- 
mission of thin film on a substrate 



In this section, we calculate the total transmission of a thin film of thickness d 
and refractive index n—ik on a transparent substrate with index s. We show that 
our result agrees with the expression given in Eq. (4) of Ref. [T] . The derivation 
is for the case where the film is weakly absorbing, so there is attenuation of light 
passing through the film (given by the attenuation factor, x), but the imaginary 
part of the refractive index, fc, is neglected in the expressions for the reflection 
and transmission coefficients at the interfaces. 




Figure SI: Definitions of the amplitude refiection and transmission coefficients 
for the film calculation. 

We first evaluate Tf and Rf for the thin film. They are obtained, in the 
usual way, using a similar approach to our derivation of Eqs.(l) and (2), except 
adding amplitudes because the reflections are coherent. On each double pass 
through the fllm, there is a phase change 4> and an amplitude reduction, a;, 
which are given by 

Airnd ( ATrkd\ 
<p = — — X = exp I — I = exp (—ad). (SI) 

The transmittance of the fllm is then 

" s|l-r2ficce**|2' ^^^^ 

where the amplitude reflection and transmission coefficients at the interfaces, 
ti, ti, t2, t2, Ti, fi, r2 and f2, are defined in Fig. (SI). Using the the Fresnel 
formulae, they are 

_ 2n ~ _ 2 _ 2s ~ _ 2n 

ti — — r^, ti — — r^, h — — \ — , h — — , — , 
n+1 n+1 n+s n+s 

n-l , s-n , 

n = — — =-ri, "^2 = — — = -7-2- (S3) 
n + 1 s + n 
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Substituting these expressions into Eq.(S2), we obtain 
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where the denominator is 

A = (n + lf{n + sf - 2x{n^ - l){n'^ - s^) cos(^ + x'^{n - lf{n - sf. 
The reflectance of the film from the substrate side is 
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where we have used the identity 
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Substituting from Eq.(S3), we get 
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where A is again given by Eq.(S5) and 

N = {n+ if (n-sf - 2x{n^ - l){n'^ - s^) cos (j) + x'^ {n - if {n + sf . (S9) 

The back transmittance and reflectance are obtained from the Fresnel for- 
mulae for the substrate-air interface: 
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We are now in position to evaluate total transmittance using Eq.(l), with a = 0. 
We have 
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Substituting A and N, the denominator of this expression beeomes 

(s + lf[{n + lf{n + sf - 2x[n^ - X){n^ - s^) eos + a;^(n - \f{n-sf\ 

- (s - \ f{{n + \f[n - sf ~ 2x{n^ - l){n^ - s^) cos0 + x'^{n - lf{n + s)^]. 

(S12) 

We evaluate the powers of x separately. First, for x": 

(s + + l)^(n + s)^ - (s - + - s)^ = 4s(n + 1)^(71 + s^) 

= 4sB, (S13) 

for x^: 

-2x{n^ - l){n^ - 5^)008 [(s + 1)2 - (s - 1)^] = -8sx(n2 _ i)(n^ _ s2)cos</> 

= -4sCxcos</) (S14) 

and for x"^: 

X^{n ~ 1)2 [(s + 1)2(„ - s)2 _ (s _ l)2(^ + 5)2] = 4,^2 _ ^)3(„ _ ^2) 

= A,sDx^, (S15) 

where the coefficients S, C and D are the same as in Ref.[T]. Putting these 
parts back together, and defining A = IGn^s, Eq.(Sll) becomes 

^ ^ B - Cx cos (j) + Dx^' ^^^^^ 

which is Eq.(4) of Ref.d]. 
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